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Abstract. Hard processes are a well calibrated probe to study heavy- ion collisions. 
However, the information to be gained from the nuclear suppression factor Raa is 
limited, hene one has to study more differential observables to do medium tomography. 
The angular correlations of hadrons associated with a hard trigger appear suitable as 
they show a rich pattern when going from low px to high px- Of prime interest is the 
^\ . fate of away side partons with an in- medium pathlength 0(several fm). At high px 

the correlations become dominated by the punchtrough of the away side parton with 
■ subsequent fragmentation. We discuss what information about the medium density 

, can be gained from the data. 
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1. Introduction 



Considerable effort is made in the study of the energy loss of hard partons in the hot 
and dense medium created in ultrarelativistic heavy-ion collisions. In particular medium 
tomography (i.e. gaining information about the medium density evolution) by studying 
the nuclear suppression factor Raa, the observed transverse momentum spectrum of 
hard hadrons divided by the scaled expectation from p-p collisions pQ. However, as 
argued in [2], Raa exhibits very limited sensitivity to the energy loss mechanism or 
properties of the medium beyond the fact that the quenching is substantial. 

If one makes the ansatz that the local quenching power q of the medium can 
be written as a function of the local energy density e as [3] q(r,r) s ,r, (ft) — K ■ 2 • 
[e(r, i] s , r, 0)] 3//4 where r is evolution proper time, r/ s spacetime rapidity and r radius of 
the medium, adjusting the scale parameter K parametrizes both the coupling strength 
of the parton to the medium and uncertainties in the determination of the medium 
evolution (see [I]). In [5] it was shown that there is about 50% uncertainty in the 
extraction of the medium quenching power (given by K) from Raa even for model 
evolutions which describe bulk soft matter properties; a 2+ld hydrodynamical evolution 
[6], a 3+ld hydrodynamical evolution [7] and a parametrized evolution [8]. Greater 
sensitivity to the medium density evolution can hence only be gained in more differential 
probes such as 7-hadron correlations [2], Raa for non-central collisions [5 J or back-to- 
back hadron correlations pE], here we focus on the latter. 
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2. The model description 

We make a Monte Carlo (MC) simulation of hard back-to-back processes in a medium in 
order to calculate the induced hard hadronic correlation pattern. The model is described 
in detail in [4]. There are three main elements: 1) the primary hard pQCD process 2) 
the description of the soft medium and 3) the energy loss from hard to soft degrees of 
freedom coupling 1) and 2). 

The energy loss for a given parton path is described in a probabilistic language. 
In order to determine the probability P(AE, E) pat f l for a hard parton with energy E 
to lose the energy AE while traversing the medium on its trajectory, we make use of a 
scaling law [9] which allows to relate the dynamical scenario to a static equivalent one 
for each trajectory £(r) by calculating 



as a function of the jet production vertex r and its angular orientation 0. We 
set q = whenever the decoupling temperature of the medium T — Tp is reached. 
Using the numerical results of [10J, we obtain P(AE) path for uj c and R = 2u^/ (qL) as a 
function of jet production vertex and the angle 4>. 

From this, the probability distribution of energy loss reflected in single hadron 
observables (such as Raa) can be computed by averaging over all possible production 
points determined by the nuclear overlap as 



The averaged energy loss probability underlying dihadron correlations involves 
a conditional probability of observing a high pj- trigger hadron, hence it is more 
complicated and must be computed in the full MC simulation. We do the simulation 
generating hard pQCD back-to-back events inside the evolving bulk matter model, using 
Eq. ((Tj) to determine P(AE) path for each parton, followed by leading and next to leading 
fragmentation once the parton emerges from the medium [I] . In all simulations, the scale 
parameter K is adjusted such that the measured Raa is described well, thus we are only 
interested in information beyond what can be extracted from single hadron observables. 

We test six different spacetime evolution scenarios: 1) a 2+ld Bjorken 
hydrodynamics simulation ('hydrodynamics') 2) a 2+1 d Bjorken hydrodynamics 
evolution in which quenching occurs only in the quark-gluon plasma phase ('black 
core'), 3) a parametrized non- Bjorken expansion with box density 4) a parametrized 
non-Bjorken expansion with nuclear profile Ta density 5) a Bjorken expansion with box 
density and 6) a Bjorken expansion with Ta density [4]. In this way, we represent a 
selection of different longitudinal and transverse dynamics. 




and 
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(2) 
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3. Results 



We compare the computed yield of associated hadrons as a function of momentum bin 
per trigger in the range 8 < px < 15 GeV on near and away side with the STAR data 
[TT| [12] in Fig. [TJ The near side is described well and there is broad agreement across 
all different medium models. The reason is that in all cases P(AE)t aa is dominated by 
either transmission without energy loss or complete absorption, the kinematic window 
to observe shifts in parton energy is very narrow. However, the probability distribution 
of experiencing a shift AE in energy is where tomographical information resides, thus, 
all scenarios appear similar for a small difference between trigger and associate px- 
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Figure 1. Yield per trigger on the near side (left) and away side (right) of hadrons 
in the 4-6 GeV and 6+ GeV momentum bin associated with a trigger with 8 < pr < 
15 GeV for different evolution models compared with STAR data [HUE]- Individual 
data points are spread artificially along the x axis for clarity. 



Since the away side yield is determined by a different (conditional) probability 
P(AE) trigger the models exhibit more pronounced differences. All models fail in the 
momentum window between 4 and 6 GeV to describe the data. This is presumably due 
to the fact that pQCD and fragmentation is not the dominant mechanism of hadron 
production in this momentum regime. Within errors, all scenarios except the 'black 
core' one agree with the data above 6 GeV associate momentum. 

The black core scenario is characterized by a very high quenching power (48 
GeV 2 /fm at r = 1 fm/c) in the center and a large halo out of which partons can escape 
unhindered. On the other hand, the hydrodynamics scenario shows less quenching power 
(11 GeV 2 /fm at r = 1 fm/c) but has a larger extension. We show the origin of dihadrons 
in each model in Fig. [2j It is evident that dihadron production in the black core scenario 
is peripheral while in the hydrodynamics scenario it occurs in the whole medium. From 
the fact that only one of these scenarios ('hydrodynamics') describes the data, we can 
infer that the produced matter is not characterized by an opaque core and a dilute halo 
and that dihadron emission is not peripheral but probes indeed the central region of the 
medium. While this is some promising start, a larger lever arm in trigger p? would be 
needed to be sensitive to shifts in parton energy instead of just absorption and to get 
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Figure 2. Probability density for finding a vertex at (x, y) leading to a triggered event 
with 8 < pr < 15 GeV and an away side hadron with 4 < px < 6 GeV for different 
spacetime evolution scenarios. The near side hadron propagates to the — x direction. 



more tomographical information [4j. Recently, elastic energy loss has been suggested 
as an additional ingredient to describe the nuclear suppression [13]. To the degree that 
the magnitude of this contribution is parametrically given by AE ~ J °° d£q(£), these 
results place tight limits on the importance of collisional energy loss: Since the integral 
has virtually no contribution from later times (long paths), elastic energy loss cannot 
resolve spatial structures larger than ~ 2 — 3 fm whereas the observed difference between 
near and away side yield relies on the fact that pathlength differences of the order of 5 
fm or more are resolved. 
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